HEAVY QUARK PRODUCTION IN pp COLLISIONS 
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Abstract 



A systematic study of the inclusive single heavy quark and heavy-quark 
pair production cross sections in pp collisions is presented for RHIC 
and LHC energies. We compare with existing data when possible. The 
dependence of the rates on the renormalization and factorization scales 
is discussed. Predictions of the cross sections are given for two different 
sets of parton distribution functions. 
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INTRODUCTION 

Charm and bottom quark production from initial nucleon-nucleon collisions will 
be copious at the RHIC and LHC colliders. Heavy quark decay into leptons will 
represent a significant background to dilepton production || in heavy ion collisions. 
A quantitative knowledge of the production cross section in pp collisions is a pre- 
requisite for the detection of collective effects, such as heavy quark production by 
rescattering and in the quark-gluon plasma, which appears as a deviation from the 
simple superposition of hadronic collisions. 

The lowest order (Born) calculations of the total cross section predict the correct 
energy dependence but differ from the experimental measurements by a U K factor" 
of 2-3. While the single-inclusive distributions as well as the mass and rapidity 
distributions of QQ pairs are also well described to within a K factor by the Born 
cross section, the pt and azimuthal double-differential distributions are not calculable 
at the Born level since the QQ pair is always produced back-to-back in lowest order. 
For this reason, a next-to-leading order (NLO) calculation is needed. The calculations 
we present here are done using a Monte Carlo program developed by Nason and 
collaborators 0, [3], |j]. Similar work on the total cross section and the single inclusive 
distributions was done by Smith, van Neerven, and collaborators ||. 

In this calculation, in addition to the uncertainties in the parton distribution 
functions, uncertainties arise from the heavy quark mass and the renormalization 
and factorization scale parameters. At collider energies, the calculations become 
more uncertain due to the lightness of the heavy quark compared to the center of 
mass energy, ttiq/ 'y/s <C 1. We first discuss the Born calculation in some detail 
and then outline the NLO calculation with its additional uncertainties. We use the 
available data on c^f(s) to fix the charm quark mass and the scale parameters. 
The resulting parameter set provides a point from which to extrapolate to heavy- 
ion collider energies. We then compare with single-inclusive and double-differential 
distributions from charm and bottom data when available. We present estimates of 
heavy quark production cross sections in proton-proton collisions at RHIC (y/s = 200 
and 500 GeV) and LHC (y/s = 5.5 TeV and 14 TeV), according to our present 
theoretical knowledge. We provide both the Born and NLO results for the total 
QQ production cross section, single inclusive y and pr distributions, and double 
differential M, 0, y and px distributions. 



HEAVY QUARK PRODUCTION IN PERTURBATIVE QCD 

The most general expression for the double differential cross section for QQ pair 
production from the collision of hadrons A and B is 

Here A and B represent the initial hadrons and i, j are the interacting partons, 
and the functions Fi are the number densities of gluons, light quarks and antiquarks 



evaluated at momentum fraction x and factorization scale fj,p. The short-distance 
cross section, a^, is calculable as a perturbation series in q> s (hr) where the strong 
coupling constant is evaluated at the renormalization scale //r. Both scales are of the 
order of the heavy quark mass. At leading order, /j, F = ft R = \i where /i = 2m c is 
commonly used. The scale dependence will be discussed in more detail below. 



Leading Order 

At leading order, O(al), QQ production proceeds by two basic processes, 

q + q -»■ Q + Q (2) 

g + g -► Q + Q . (3) 

The invariant cross section for the process A + B — > H + H where the QQ pair has 
fragmented into hadrons H(Qq) and H(Qq) can be written as 
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where y§, the parton-parton center of mass energy, is related to y/s, the hadron- 
hadron center of mass energy, by s = X1X2S. The intrinsic transverse momenta of the 
incoming partons have been neglected. The sum of the leading order subprocess cross 
sections convoluted with the parton number densities is contained in C(x\,X2) where 

C(x 1 ,x 2 ) = £[i?(si)F*(s 2 ) + if (xi)if (* 2 )]^ + F*{ Xl )Ff(x a )^ • (5) 

Only light quark flavors, those with m < ttiq, are included in the sum over q. The 
dependence on the scale \i F has been suppressed here. 

Fragmentation affects the charmed hadron distributions, not the total cc produc- 
tion cross section. The fragmentation functions, Dh/q{z), describe the hadronization 
of the heavy quarks where z = \ph\/\Pq\ is the fraction of the heavy quark momen- 
tum carried by the final-state hadron. The D meson x F distribution is harder than 
the calculated charmed quark distribution in hadron-hadron interactions. Including 
a fragmentation function that describes D production in e + e~ annihilation softens 
the distribution due to energy lost to light qq pair production ||. Event genera- 
tors such as PYTHIA [0, based on the Lund string fragmentation model, harden 
the D distribution. In PYTHIA, the charmed quark is always at the endpoint of a 
string which pulls the charmed quark in the direction of a beam remnant so that the 
charmed hadron can be produced at a larger momentum than the charmed quark. 
Correlations of the produced charmed hadron with the projectile valence quarks, not 
predicted by perturbative QCD, have been measured. Several possible explanations 



have been suggested, see i.e., |6|, [7|, §]. This interesting high ii? regime will not be 
measurable at the RHIC and LHC colliders since the center of mass energy is high 
and the rapidity coverage is mostly confined to the central region. (The PHENIX 
muon spectrometer at RHIC will have a larger rapidity coverage, 1.5 < y < 2.5 ||, 
but these effects will probably be out of reach at the maximum collider energy.) 

If we ignore fragmentation effects for the moment, after taking four- momentum 
conservation into account, we are left with 



x 1 x 2 C{x ll x 2 ) , (6) 



dp 2 T dy Q dy-Q 



where X\ and x 2 are 



Xl = ^(e^ + e^) , (7) 

X2 = H^( e -yQ +e -^) 

and m,Q = J frig + p 2 -,. At yq = y-Q = 0, x\ = x 2 . The target fractions, x 2 , decrease 
with rapidity while the projectile fractions, X\, increase. The subprocess cross sections 
for QQ production by qq annihilation and gg fusion to order 0(a 2 ), expressed as a 



function of uiq, |/q, and y-g are fT0| 

doqq = na 2 s cosh(y Q - y$) + m 2 Q /m 2 Q 
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da gg _ ira 2 8 cosh(y Q - y^) - 1 f lf N , 2m 2 Q 2m, 



di 96m| (1 + cosh(y - y^))' 



{0S h(y Q -yQ) + ^--^\ . (9) 



Next-to-Leading Order 

We now discuss the NLO, 0(a?), corrections to the QQ production cross sec- 
tion. At next-to-leading order, in addition to virtual corrections to these diagrams, 
production by 

(10) 

(11) 
(12) 

must also be included. The last process, quark-gluon scattering, is not present at 
leading order. The quark-gluon graphs can be interpreted at the Born level as the 
scattering of a heavy quark excited from the nucleon sea with a light quark or gluon 



q + q - 


-> Q + Q + g 


9 + 9 ~ 


-> Q + Q + g 


Q(Q) +9 ~ 


-> Q + Q + (q)q 



and are referred to as flavor excitation [0. The total short distance cross section d^ 



for a given production process can be expressed generally as 

a ij (s,m Q ,p R ) = s - 2 R fij{p,li 2 R /m 2 Q ) , (13) 

m Q 

where p = irrig/s. The function fy can be expanded perturbatively as 

UP, A/™%) = 4(p) + ^^ [flM + J),{P) HA/™%)] + 0(a 2 s ) . (14) 

The leading order part of the cross section is in the function /P. In this case, /°, = 
$9 = fgq = fg\ = ^' Only f® and f^ contribute and can be computed from the t 
integration of the cross sections given in (8) and (9). The physical cross section should 
be independent of the renormalization scale: the dependence in eq. (14) introduces an 
unphysical parameter in the calculation. If the perturbative expansion is sufficient, 
i.e. if further higher-order corrections are small, at some value of p the physical 
0(a™ ) and (D(a™ +1 ) cross sections should be equafj. If the p dependence is strong, the 
perturbative expansion is untrustworthy and the predictive power of the calculation 
is weak JTOf- The rather large difference between the heavy-quark Born and NLO 
cross sections suggests that further higher-order corrections are needed, particularly 
for charm and bottom quarks which are rather "light" when yfs is large. Usually 
the renormalization scale in a^ and the factorization scale in the parton distribution 
functions are chosen to be equal. We follow this prescription in our calculations. 

We have used two sets of recent parton distribution functions^], GRV HO [I2| and 
MRS D-' |R|. The first begins with a low scale, Qqgrv = 0-3 GeV 2 , and valence- 
like parton distributions, therefore evolving very quickly with Q 2 . The second, with 
Qomrs = 5 GeV 2 , has sea quark and gluon distributions that grow as ~ x -1 / 2 when 
x — > 0. Both are compatible with the recent deep- inelastic scattering data from 
HERA |I3] . We also include estimates of the total cross section using the MRS DO' [13| 



distributions. This set assumes a constant value for the sea and gluon distributions at 
Qo mrs as x — >■ and lies below the HERA data. The GRV distributions assume u = d, 
a symmetric light quark sea, and xs(x, Qogrv) = 0' increasing to give 2{x)s/((x)u + 
(z) 5 ) ~ 0.53 at Q 2 = 10 GeV 2 [0. The MRS D sets allow u <_ d to account 
for measurements of the Gottfried sum rule and assume s = (u + d)/A at Qomrs 
[ |T3]j . Thus the MRS distributions, arising from a global fit, provide a somewhat 
better description of the deep-inelastic scattering data for x > 0.01 than the GRV 



distributions |12l 13 



Since we compare two extreme cases for the nucleon parton distributions as x — > 0, 
MRS D-' and GRV HO on one hand and MRS DO' on the other, our results may be 
thought of as providing an upper and lower bound to the QQ cross section at heavy- 
ion collider energies for fixed mass and scale. However, little data exist on the gluon 

■tThe order of the expansion is represented by n. For QQ production, n > 2. A calculation to 
order 0(a") introduces corrections at the order 0(a" +1 ). 

''All available parton distribution functions are contained in the package PDFLIB |0|, available 
in the CERN library routines. 



distribution function at low x so that it is poorly known, particularly in the x region 
accessible at RHIC and LHC, x ~ 10~ 2 and 10 -4 around y = 0, respectively. The low 
x behavior has a significant effect on the shape of the gluon distribution at moderate 
values of x in the energy range of Fig. 1 . Steeply rising gluon distributions at low x 
are compensated for by a corresponding depletion at moderate x. 

Heavy quark production by gluon fusion dominates the pp —>■ QQX production 
cross section in the central region. Thus we show the shape of the gluon distributions 
of the three parton distribution sets are shown in Fig. 1(a) over the x range of the 
previous pp data, 0.01 < x < 1. To facilitate comparison, all three are shown at 
// = 2.4 GeV. The solid curve is the GRV HO distribution, the dashed, MRS DO', 
and the dot-dashed, MRS D-'. The GRV distribution at /i = 1.2 GeV is also shown 
to demonstrate the effect of the Q 2 evolution. Since it has a smaller initial scale, the 
evolution with \x is quite fast. The DO' distribution can be seen to turn over and 
begin to flatten as x decreases. However, for much of the range, it is above the D-' 
distribution, reflected in a larger a 1 ^, as shown in Fig. 3. All three sets, evaluated in 
the MS scheme, have a similar value of Aqcd- in Fig. 1(b), we show the running of 
the two loop value of a s , 



®s{H,f) 
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where 6/ = (33-2/)/12tt, b' f = (153-19/)/(2tt(33-2/)), / is the number of flavors, 
and Af is the value of Aqcd appropriate for the number of flavors. In the calculation, 
the number of flavors depends on the chosen quark mass. For charm, / = 3, and 
for beauty, / = 4. At [A — rriQ, a s (rriQ, f) = a s (rriQ, f + 1). The running of a s is 
visible in the renormalization scale dependence, shown in Fig. 2(e). For the NLO 
QQ production program, Aj is chosen by mq. Note that A 3 > A 4 > A 5 . Additional 
uncertainties may arise because the threshold ttiq for a given parton distribution set 
can differ from our fitted vhq. 

While it is often possible to use a general prescription like the principle of minimal 
sensitivity (PMS) |0| to find values of /i# and \ip where the scale sensitivity is a 
minimum, the heavy quark production cross section is very sensitive to changes in \i. 
In Fig. 2 we show the variation of the cc and bb production cross sections at RHIC 
(a), (c) and LHC (b), (d) ion energies. The MRS distributions exhibit an artificial 
stability for low \x because for // < 2m c ps Qo,mrs> the factorization scale is fixed at 
Qo.mrs an d only ji R varies. We use the GRV HO parton distribution functions so that 
we can show the uncertainty with \i = ji R = /x F at lower values of \i since \ip is not 
fixed until \i F ~ 0.4m c k, Qo,grv- When fi/m c w 0.2, the cross section diverges since 
(fjb/m c )/ Aqcd ~ 1. In any case, such small scales below 1 GeV, are excluded because 
a perturbative calculation is no longer assumed to be valid. As fi/m c increases, the 
cross section becomes more stable. The behavior we find is similar for RHIC and 
LHC energies. The bb cross section shows a smaller variation with fi, particularly 
at y/s = 200 GeV. The variation resembles the running of a s shown in Fig. 1(b). 
Indeed, this running is a major source of instability in the NLO QQ cross sections. 
However, at yfs = 5.5 TeV the variation with /i at the Born level increases since 



the cross section becomes more uncertain as tuq/ ' yfs decreases. The NLO results 
show less variation at this energy. There is no value of /i where the Born and the 
NLO calculations are equal, suggesting that higher-order corrections are needed for 

m Q/Vs < 1. 

We show the change of the cc cross section at yfs = 200 GeV induced by fixing 
fin = 2niQ and changing \x F in Fig. 2(e) and fixing fip = 2tjiq and varying fi R in 
Fig. 2(f). The running of the coupling constant is clearly shown in 2(e). In 2(f), the 
increase with ftp arises because at values of \ip near Qo.grv and low x, the sea quark 
and gluon distributions show a valence-like behavior, decreasing as x — > 0, an effect 
special to the GRV distributions |TJ]. The results are quite different for the MRS 
distributions, especially for the equivalent of Fig. 2(f). There is not much change in 
the cross section with /j,p, particularly at the Born level, since the parton distribution 
functions do not change below Qo,mrs- 



CALCULATIONS OF a^ 

Previous comparisons of the total charm production cross sections with calcula- 
tions [[[(J at leading order suggested that a constant K factor of ~ 2 was needed to 
reconcile the calculations with data when using m c = 1.5 GeV, but not when m c = 1.2 
GeV was chosen. Initial NLO calculations seemed to suggest that the K factor was 



no longer needed with m c = 1.5 GeV ||17||. However, this result is very dependent 



upon the chosen scale parameters and the parton distribution functions, particularly 
the shape of the gluon distribution. 



Comparison With Current Data 



on 



We compare our NLO calculations with the available data [JIB, [19], |20|, |21], |22 
the total cc production cross section from pp and pA interactions in Fig. 3. When 
a nuclear target has been used, the cross section per nucleon is given, assuming 
an A a dependence with a = 1, supported by recent experimental studies of the 
A dependence |E3|. We assume that we can compare the cc production cross sec- 
tion directly with charmed hadron measurements. Often single charmed mesons, 
denoted D/D to include all charge states, in the region xp > are measured. The 
cc production cross section is symmetric around xp = in pp interactions so that 
a^ = 2a C c(xp > 0). While the question of how the cc pair hadronizes into DD, DA C , 
A C D, A C A C , etc. remains open, some assumptions must be made about how much of 
cr^f is missing since not all channels are measured. If all single D mesons are assumed 
to originate from DD pairs, ignoring associated A C D production, then by definition 
cr(DD) = a(D/D)/2. Thus the single D cross section for xp > is equal to the DD 
pair cross section over all xp. However, the contribution to the cc total cross section 
from D s and A c production has been estimated to be a(D s )/a(D° + D + ) ~ 0.2 and 
cr(A c )/cr(/}° + D + ) ~ 0.3. Thus to obtain the total cc cross section from cr(DD), 



a(DD) should be multiplied by ~ 1.5 ^2M. This is done in our data comparison. The 
data exist in the range 19 < \fs < 63 GeV, mostly from fixed target experiments. 
Below the ISR energies, y'i = 53-63 GeV, the total cross section is primarily inferred 
from single D or DD measurements. At the ISR, the pair production cross section is 
obtained from lepton measurements, either e/x and electron pair coincidence measure- 
ments or a lepton trigger in coincidence with a reconstructed D or A c . Rather large 
cc cross sections were inferred from the latter analyses due to the assumed shape of 
the production cross sections: flat distributions in xp for the A c and (1 — xp) 3 for the 
D. The ISR results must thus be taken with some care. 

Modern parton distributions with Aqcd fixed by fits to data cannot explain the 
energy dependence of the total cross section in the measured energy range when using 
m c = 1.5 GeV and \ip = Hr = m c . Since m 2 < Ql MRS for the MRS distributions and 
the scale must be chosen so that ji 2 > Ql MRS for the calculations to make sense, we 
take fj, = 2m c and vary m c for these distributions. We find reasonable agreement for 
m c = 1.2 GeV for the D-' and DO' distributions. The results are shown in the solid and 
dashed curves in Fig. 3 respectively. Since the GRV HO distributions have a much 
lower initial scale, \x can be fixed to the quark mass. The dot-dashed curve is the GRV 
HO distribution with m c = 1.3 GeV and \x = m c . All three curves give an equivalent 
description of the data. Our "fits" to the low energy data are to provide a reasonable 
point from which to extrapolate to higher energies. It is important to remember that 
significant uncertainties still exist which could change our estimates considerably 
when accounted for. These relatively low values of m c effectively provide an upper 
bound on the charm production cross section at high energies. For comparison, we 
also show the cross section with the GRV distributions and \x = m c = 1.5 GeV in the 
dotted curve. It lies a factor of 2-3 below the other calculations. The smaller value 
of m c is needed for the MRS distributions even with the larger scale because parton 
distribution functions at lower values of Q 2 would decrease at low x, as demonstrated 



by the GRV distributions p2[ . Note that such small choices of m c suggests that the 
bulk of the total cross section comes from invariant masses less than Imp,. In a recent 
work [23], the total cross section data was found to be in agreement with m c = 1.5 
GeV with some essential caveats: the factorization scale was fixed at fip = 2m c while 
fin was allowed to vary and an older set of parton distribution functions with a range 
of fits with a different value of Aqcd for each was used. Decreasing /i# with respect to 
fip and increasing Aqcd both result in a significantly larger cross section for a given 
m c . We choose here to use the most up-to-date parton distribution functions and to 
keep fip = fiR, facilitating a more direct extrapolation from the current data to the 
future collider results. 

Since data on cc and bb production by pion beams are also available at fixed 
target energies, in Fig. 4 we show this data with the same parton distributions where 

is 



m c and /i are fixed by the comparison in Fig. 3. The cc data [18|, 25| |26|, £7], |28 



based on the xp > single D cross section. However, the n~N xp distribution is 
asymmetric, o jo{xp > 0) ~ 1.6 so that a(DD) is obtained by dividing by 2 to get the 
pair cross section and then multiplying by 1.6 to account for the partial xp coverage. 
The bb data, taken to be over all xf, are generally obtained from multi-muon studies 



[29 , |30| , [31 , [32| . The data, especially for bb production, are not as extensive and have 



rather poor statistics. Again, some of the data is from a nuclear target. When a 
nuclear target has been used, the cross section per nucleon is given, assuming an A 1 
dependence. 

The GRV HO pion distributions |33| are based on their proton set so that the 
two distributions are compatible. In Fig. 4(a), the charm production cross section is 
calculated using the GRV proton and pion distributions. The solid curve shows the 
result with a nucleon target, the averaged distributions for proton and neutron, while 
the dashed curve is the result for a proton target alone. The results are consistent at 
yfs = 30 GeV; at lower energies, the cross section on a proton target is slightly larger 
than on a nucleon target. The calculations using the MRS distributions do not have 
the same consistency as those with GRV because their pion distribution functions, 



SMRS PI and P2 |3"^{] , are based on an older set of proton distributions than the 
current MRS distributions used here. The SMRS distributions use A 4 = 190 MeV 
while the MRS distributions have fixed A4 = 230 MeV. In the calculations, we fix A 4 
to the current MRS value. The dot-dashed curve shows the MRS D-' distributions 
with the SMRS P2 pion distributions while the dotted curve is with the PI set. Both 
are for a proton target. The PI set has a steeper gluon distribution than P2. The 
two calculations begin to diverge as a/s increases since the gluon fusion contribution 
is becoming dominant. At low \/s, valence quark annihilation is important for n~p 
interactions. Although the calculations and data are not in exact agreement, they 
are close enough to assume that the same parameters are reasonable for both pion 
and proton projectiles. The comparison to the bb production cross section is given in 
Fig. 4(b). The data is very sparse. We use m b = 4.75 GeV and \x = mb for both sets 
of parton distributions. The solid curve is the GRV distribution, the dashed is the 
MRS D-' and SMRS PI result. The agreement is not unreasonable given the quality 
of the data on the one hand and the theoretical uncertainties on the other. 



Extrapolation To RHIC And LHC Energies 

The total cc and bb cross sections at the top ISR energy, a/s = 63 GeV, and 
the proton and ion beam energies at RHIC and LHC are given in Tables 1 and 2 
respectively. Both the Born and NLO cross sections are given. The theoretical K 
factor, fng /^, tends to increase with energy and is rather large. There is no a 
priori reason why it should remain constant, rather the increase at collider energies 
would suggest that the perturbative expansion is becoming less reliable, as discussed 
below. Note that even though the MRS D-' and GRV HO distributions give an equally 
valid description of the data at ISR energies and below, they differ at higher energies, 
partly from the difference in m c and partly because of our scale difference. The MRS 
D-' distributions evolve faster since fi = 2m c rather than /i = m c due to their chosen 
initial scale Qo,mrs 5 resulting in a larger predicted cross section. Less difference is 
seen between the GRV and MRS D-' distributions for the bb cross section since the nib 
and \i are used for both. Note that for bb production at 14 TeV, the results differ by 
30% while the MRS D-' NLO cc result is three times larger than the GRV HO result 



at the same energy. The DO' distributions give smaller cross sections at LHC energies 
due to the different initial behavior at x — > 0. We illustrate this effect using the Born 
contribution to the production cross section at fixed M and y = 0, approximated as 



da a 2 r ,- 2 



lj/=o ■- 



dMdy y Ms 



F 9 (M/V~s) (16) 



since gluon fusion is the dominant contribution to the Born cross section, x = Mfy/s 
at y — 0, and at fixed M, a gg is proportional to (a 2 /M 2 )F 2 . The gluon distribution 
at low x and fi = Qo may be approximated as F g (x) = f(x)/x 1+s . For a constant 
behavior at low x, such as in the MRS DO' distribution, 5 = and the cross section is 
independent of ^Js. At the other extreme, the MRS D-' distribution assumes S = 0.5 
at Qq so that the cross section grows as s 5 ~ */s. 



SINGLE AND DOUBLE DIFFERENTIAL DISTRIBUTIONS 

We now compare the NLO calculations with data on Q and QQ distributions. 
In the presentation of the single inclusive and double differential distributions, we 
follow the prescription of Nason and collaborators |3|, f| and take us = nffiQ for 
the single and fxjj = nJrriQ + (p 2 -. +Pt_ )/2 for the double differential distributions. 
When using MRS distributions for charm production, n = 2. For all other cases, 
n = 1. A word of caution is necessary when looking at our predictions for QQ pair 
distributions. It is difficult to properly regularize the soft and collinear divergences to 
obtain a finite cross section over all phase space. Soft divergences cancel between real 
and virtual corrections when properly regularized. The collinear divergences need to 
be regularized and subtracted. For single inclusive heavy quark production, this is 
possible because the integration over the partonic recoil variables can be performed 
analytically and the singularities isolated. In exclusive QQ pair production, the 
cancellation is performed within the numerical integration. The price paid for this 
is often a negative cross section near the phase space boundaries, particularly when 
Pt — > for the pair and <fi — > it where is the difference in the azimuthal angle 
between the heavy quark and antiquark in the plane transverse to the beam axis. 
A positive differential cross section for p? — > can only be obtained by resumming 
the full series of leading Sudakov logarithms corresponding to an arbitrary number of 
soft gluons. This has not been done in the case of heavy quark production [||. Thus 
when mq/^/s <^. 1, fluctuations in the cross section due to incomplete numerical 
cancellations can become very large, resulting in negative components in the mass 
and rapidity distributions. We have minimized the fluctuations by maximizing the 
event sampling at low px and increasing the number of iterations |35[. 



Comparison To Current Data 



First, we compare with the 800 GeV fixed target data of the LEBC-MPS col- 
laboration [19| in Fig. 5. They measured the xf and p T distributions of single D 



production. The total cross section, a(D/D) = 48 ± 11 /ib, corresponds to a DD 
production cross section of 24±8 /ib. The solid curves are the MRS D-' results, the 
dashed, the GRV HO calculations. Data on correlated DD production is also avail- 



able at 800 GeV, from pEmulsion studies |p6fl . The event sample is rather small, only 
35 correlated pairs. We compare the mass and p T of the pair and the azimuthal dif- 
ference between the pair in Fig. 6 with the calculated NLO distributions. Again the 
solid curve is MRS D-', the dashed, GRV HO. The Born invariant mass distribution, 
given by the dashed curve, is parallel to the NLO results shown in the solid curve. 

The pp data from UA 1, y/s = 630 GeV, and CDF, y/s = 1.8 TeV, include 
single b quark px distributions. The measurements are taken in the central region 
(\y\ < 1.5 for UA 1 and \y\ < 1 for CDF) and are integrated over px above each 
Pr,min- The comparisons with the NLO calculations are given in Fig. 7(a) for UA 1 
[[37] and Fig. 7(b) for CDF and DO [38, 3J5 . Reasonable agreement is found for both 



GRV HO and MRS D-' for UA 1 with [is — \f m b + Pt- However, the results from 
this same scale choice lie somewhat below the early CDF data where data on J/ip 
production was used to determine the B production cross section^. As reported in 
Ref. |40| , the scale \x = /is/4 was needed for good agreement with the magnitude of 



the data when the older MRS DO distributions were used. More recent data using 
direct measurement of inclusive b — > J/ip and b — > tp' decays has shown that the 
previous results overestimated ip production from b decays fl58"| . Better agreement 
with theory is now found for \x = [is, as shown in Fig. 7(b). Again the GRV HO 
and MRS D-' distributions look similar, differing primarily for pr,min < 10 GeV. This 
difference is increased for the lower scale choice where /is/4 < Qo,mrs for px,mm < 7.5 
GeV, cutting off the evolution of the MRS distributions below this px,min- The GRV 
calculations evolve over all PT,min since /is/4 > Qo,grv> hence the larger difference. 



Extrapolation To RHIC And LHC Energies 

We now show the predicted heavy quark distributions for RHIC (\/s = 200 and 500 
GeV) and LHC (y/s = 5.5 and 14 TeV) using the MRS D-' and GRV HO distributions. 
The results are shown in Figs. 8-23. We use the same scales on the ?/-axes for both sets 
of parton distributions as much as possible to facilitate comparison. In each figure 
we show the single quark px (a) and y (b) distributions and the px (c), rapidity (d), 
invariant mass (e), and azimuthal angle (f) distributions of the QQ pair. The Born 
(LO) results are also given in (b), (d), and (e). All the distributions have been divided 
by the corresponding bin width. The single and pair px distributions are also given 
with the rapidity cuts y < |1| at the LHC and y < |0.35| at RHIC, corresponding 
to the planned acceptances of ALICE |nj and the PHENIX central detector ||. 



These px distributions are also divided by the width of the rapidity interval. In 



^The inclusive decay, B — > J/tpX, has a 1% branching ratio (BR) while the channel B — * J/ipK 
has an 0.1% branching ratio. 



Tables 3-10 we give the ^/-integrated single p T NLO and Born distributions, the pair 
p T distributions with the cut on rapidity, and the NLO and Born invariant mass 
distributions for c and b production at each energy with the MRS D-' partons. Note 
that all distributions have a 2 GeV bin width and that neither it nor the rapidity bin 
width has been removed in the tables. The statistical uncertainties are less than 1% 
at low p T and M, increasing to 5-6% in the tails. The uncertainty increases slightly 
with energy. 

The development of a rapidity plateau can be seen in both the single and pair 
rapidity distributions as the energy increases. This plateau is generally broader for 
the single quarks than the pair since the pair mass enters into the estimate of the 
maximum pair rapidity while the smaller quark transverse mass gives the maximum 
single quark rapidity. The plateau is broader for the MRS D-' parton distributions. 
In the charm rapidity distributions with the MRS D-' partons at 14 TeV, the plateau 
edge is artificial. The set has a minimum x of 10 -5 , reached at y ~ 2.8 for a single 
quark and a somewhat larger y for the pair. The GRV HO distributions have a 
minimum x = 10 -6 , corresponding to y ~ 4.5, off the scale of our graphs. The 
average quark and pair px increases with energy. For charmed quarks, (p T ) is larger 
for the pair than for a single quark. The opposite result is seen for b quarks. The 
GRV distributions result in larger (p T ) than the MRS distributions. Near px — > 0, 
the MRS parton distributions show a steeper slope than the GRV distributions. As 
Pt increases, the slopes become somewhat similar at RHIC energies. 

In general, the LO mass and rapidity distributions are nearly equivalent to the 
NLO results scaled by a theoretical K factor independent of M and y. At LHC en- 
ergies, the expansion parameter becomes a s log(s/mg), of order 1 for mq/^fs <C 1, 
spoiling the convergence of the perturbative expansion |M| . This causes our predic- 
tions to be less reliable at these energies. Note that using n$ for the single inclusive 
distributions and /i£> for the double differential distributions leads to somewhat dif- 
ferent values of the integrated NLO cross sections than given in Tables 1 and 2, 
calculated with /i = nrnq, since the correction terms grow with fi. The effect is rela- 
tively small for the Born results since the faster evolution of the parton distribution 
functions is partly compensated by the decrease of a s with increasing \x. 

We also compare to the leading order charm distributions obtained from HIJING 
H for the ion collider energies, 200 GeV (Figs. 8,9) and 5.5 TeV (Figs. 12,13). 
HIJING uses the same mass and scale parameters and parton distribution functions 
as the other calculations. Although only a Born level calculation of QQ production, 
HIJING includes the effect of multiple parton showers which simulates aspects of 
higher-order production (NLO includes the effect of only one additional parton). The 
rapidity distributions, shown for y > only, closely resemble the NLO calculations. 
However, the p T distributions, taken in the rapidity interval \y\ < 2 for the single c 
quark and the pair, are softer, especially for the cc pair. (Again, the distributions are 
divided by the rapidity bin width.) The distributions are also not strongly peaked 
at low pt, as are the NLO calculations, due to initial state radiation. HIJING also 
includes fragmentation of the cc pair into hadrons. The calculated <fi distributions are 
not as sharply peaked at = it as the NLO results. Note also that the DD pair (f) 
distributions from HIJING are more isotropic than the original cc pairs. 



QQ Decays To Lepton Pairs 

Since heavy quark decays are an important contribution to the dilepton contin- 
uum, we show cc and bb decays into dileptons at RHIC and LHC for the MRS D-' 
sets. Because heavy quark decays are not incorporated into our double-differential 
calculation, the heavy quark pairs have been created from the final distributions. The 
heavy quark decays to leptons are thus calculated using a Monte Carlo program based 
on data from D decays at SLAC K3J and B decays from CLEO |44j|. The inclusive 



branching ratio for D meson decay into a lepton, averaged over charged and neutral 
.D's is BR(D°/D + — > l + X) ~ 12%. The corresponding branching ratio for B mesons 
of unspecified charge is BR(B — > l + X) ~ 10.4% [j45"| . B decays represent a special 



challenge since lepton pairs of opposite sign can be produced from the decay of a 
single B by B — ► DIX followed by D — > IX. Thus the B decays can produce dilep- 
tons from the following: a combination of leptons from a single B, two leptons from 
primary B decays, two leptons from secondary decays, and a primary lepton from 
one B and a secondary lepton from the opposite sign B. The measurement of Ref. 
|HJ ] is assumed to be for primary B decays to leptons. The NLO pair distributions 
da/dM and da/dy agree well with a K factor times the Born results. Therefore the 
correlated distributions, da/dMdy, are calculated at leading order and multiplied by 
this K-factor, while the p\ and <fi distributions, unavailable at leading order, are taken 
from the NLO results. The heavy quark pair is specified according to the correlated 
distributions from the calculated cross section. The momentum vectors of the indi- 
vidual quarks are computed in the pair rest frame, using the rapidity gap between the 
quarks. Once the quark four-momenta have been specified, the decays are calculated 
in the quark rest frame, according to the measured lepton momentum distributions, 
and then boosted back to the nucleon-nucleon center of mass, the lab frame for RHIC 
and LHC. Finally, the pair quantities, Mu, yu, and Pt,u, are computed. 

The average number of QQ pairs, Nqq, produced in a central nuclear collision 
is estimated by multiplying the cross section from Tables 1 and 2 by the nuclear 
thickness Iab(0). If ^qq < 1> onr y correlated production is important. The number 
of correlated lepton pairs can be estimated by multiplying the number of QQ pairs 
by the square of the meson, H, branching ratio to leptons: NqqBR 2 (H/H — > IX). 
However, if Nqq > 1, dilepton production from uncorrected QQ pairs should be 
accounted for as well. Then two QQ pairs are generated from the production cross 
section and the Q from one pair is decayed with the Q from the other. Thus for 
uncorrelated QQ production, the average number of lepton pairs is approximately 
N Q q(N Q q - 1)BR 2 (HH -> l±X) when N Q q > 1. If N Q q « 1, a distribution in 
Nqq must be considered to calculate the uncorrelated pairs. In the following figures, 
we show the correlated dilepton cross section in pp collisions, an = BR 2 (H/H — > 
1 ± X)<Jqq. In Fig. 27, showing uncorrelated lepton pairs from DD decays at the LHC, 
we give the uncorrelated distributions with the value of the correlated cross section 
since Nqq < 1 in pp collisions. To find the correct scale in central AB collisions, 
calculate Nqq and then multiply the lepton pair cross section by Tab (0) {Nqq — !)• 

In Figs. 24-25, we show the mass (a), rapidity (b), and pt (c) distributions for the 



lepton pairs from DD and BB pairs respectively. The average mass of the lepton pairs 
from DD decays at RHIC ion energies is {Mu) = 1.35 GeV and the average lepton 
pair pr, {p T:U ) = 0.8 GeV; from BB decays, (M«) = 3.17 GeV and {p T ,u) = 1.9 GeV. 
A like-sign subtraction should eliminate most of the uncorrelated charm production 
at RHIC. 

At LHC ion energies, the cc production cross sections are large enough for un- 
correlated charm production to be substantial and difficult to subtract in nuclear 
collisions. The average mass of the lepton pairs from correlated DD decays here is 
(Mu) = 1.46 GeV and the (pt,u) = 0.82 GeV. When the pairs are assumed to be 
uncorrelated, then {Mu) = 2.73 GeV and {pr,u) = 1 GeV. The average dilepton mass 
from uncorrelated DD pairs is larger since the rapidity gap between uncorrelated D 
and D mesons is larger on average than between correlated DD pairs. The bb cross 
section is still small enough at the LHC for uncorrelated lepton pair production from 
B meson decays to be small. However, the acceptance for these pairs will be larger 
than for charm decays since high mass lepton pairs from heavy quark decays have 
a large rapidity gap. When acceptance cuts are applied, at least one member of a 
lepton pair will have a large enough rapidity to escape undetected so that high mass 
pairs from heavy quark decays will have a strongly reduced acceptance. This reduc- 
tion will occur at larger values of Mu for BB than DD decays. From all BB decays, 
{Mu) — 3.39 GeV and (pt,u) = 2 GeV. In Figs. 26-28, we show the mass (a), rapidity 
(b), and pt (c) distributions for the dilepton pairs from correlated and uncorrelated 
DD and correlated BB pairs respectively. 



SUMMARY 

In this overview, we have attempted to use the theoretical state of the art to pre- 
dict heavy quark production in pp collisions at RHIC and LHC energies. Although 
much progress has been made in the higher-order calculations of QQ production, this 
is not meant to be the final word. Fragmentation and decay effects need to be in- 
corporated into our next-to-leading order calculations. More structure function data 
from HERA, combined with collider data on jets and prompt photons, will produce 
further refined sets of parton distribution functions. Theoretical progress may allow 
resummation at low p? or produce estimates of next-to-next-to-leading order correc- 
tions. New scale fixing techniques may result in a reduction of scale uncertainties. 
Thus, there is still room for improvement in these calculations. Though the agree- 
ment with lower energy data allows us to extrapolate these results to RHIC and 
LHC energies, major uncertainties still exist, particularly at LHC energies. However, 
given our mass and scale parameters, the GRV HO and MRS D-' parton distribution 
functions provide a rough upper and lower limit on the theoretical predictions. This 
might be useful in particular for the design of detectors at these facilities. 
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Table 1: Total cc production cross sections at collider energies. 
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Table 2: Total 66 production cross sections at collider energies. 



re- 



Production v/i = 200 GeV 



da c jdp\ 


(fib/2 GeV 2 ) 


da c - c /dp 2 T dy (/ib/2 GeV 2 ) 


da C c/dM (/ib/2 


GeV) 


p 2 T (GeV 2 ) 


NLO 


LO 


p 2 T (GeV 2 ) 


NLO 


M (GeV) 


NLO 


LO 


1 


232.5 


102.2 


1 


30.90 








3 


37.93 


15.14 


3 


3.916 


3 


172.8 


76.41 


5 


12.37 


4.589 


5 


1.548 


5 


77.05 


34.18 


7 


5.362 


1.924 


7 


0.8435 


7 


22.60 


9.611 


9 


2.774 


0.9704 


9 


0.4770 


9 


8.548 


3.429 


11 


1.589 


0.5435 


11 


0.3287 


11 


3.671 


1.427 


13 


1.003 


0.3389 


13 


0.2203 


13 


1.863 


0.6871 


15 


0.6715 


0.2206 


15 


0.1608 


15 


0.9122 


0.3438 


17 


0.4612 


0.1542 


17 


0.1277 


17 


0.5120 


0.1917 


19 


0.3291 


0.1079 


19 


0.0925 


19 


0.3154 


0.1095 


21 


0.2399 


0.0812 


21 


0.0786 


21 


0.1883 


0.0651 


23 


0.1857 


0.0602 


23 


0.0589 


23 


0.1210 


0.0415 


25 


0.1369 


0.0428 


25 


0.0478 


25 


0.0689 


0.0245 


27 


0.1088 


0.0355 


27 


0.0356 


27 


0.0520 


0.0166 


29 


0.0864 


0.0280 


29 


0.0350 


29 


0.0364 


0.0105 


31 


0.0697 


0.0225 


31 


0.0282 


31 


0.0257 


0.00785 


33 


0.0574 


0.0191 


33 


0.0206 


33 


0.0151 


0.00538 


35 


0.0478 


0.0160 


35 


0.0214 


35 


0.0111 


0.00383 


37 


0.0400 


0.0132 


37 


0.0160 


37 


0.0678 


0.00222 


39 


0.0343 


0.0111 


39 


0.0135 


39 


0.0480 


0.00198 



Table 3: 
[The rapidity-integrated p\ distribution is given for single charm (NLO and Born) 
and the p\ distribution in the range \y\ < 0.35 is given for cc pair production (NLO 
only). The tabulated results have not been corrected for the rapidity bin width. 
The rapidity-integrated pair mass distribution is also given. All distributions are at 
y/s = 200 GeV and calculated with MRS D-' parton distributions. Note the 2 GeV 
bin width for the distributions.! 
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Table 4: 
[The rapidity-integrated p\ distribution is given for single charm (NLO and Born) 
and the p\ distribution in the range \y\ < 0.35 is given for cc pair production (NLO 
only). The tabulated results have not been corrected for the rapidity bin width. 
The rapidity-integrated pair mass distribution is also given. All distributions are at 
y/s = 500 GeV and calculated with MRS D-' parton distributions. Note the 2 GeV 
bin width for the distributions.! 
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39 


4.680 


39 


3.500 


1.039 



Table 5: 
[The rapidity-integrated p\ distribution is given for single charm (NLO and Born) 
and the p\ distribution in the range \y\ < 1 is given for cc pair production (NLO 
only). The tabulated results have not been corrected for the rapidity bin width. 
The rapidity-integrated pair mass distribution is also given. All distributions are at 
y/s = 5.5 TeV and calculated with MRS D-' parton distributions. Note the 2 GeV 
bin width for the distributions.! 



cc 



Production y/s = 14 TeV 



da c /dp\ 


(fib/2 GeV 2 ) 


da c - c /dp 2 T dy (fib/2 GeV 2 ) 


da c - c /dM (nb/2 GeV) 


p 2 T (GeV 2 ) 


NLO 


LO 


p 2 T (GeV 2 ) 


NLO 


M (GeV) 


NLO 


LO 


1 


23650. 


11960. 


1 


4594. 








3 


6067. 


2473. 


3 


1129. 


3 


17250. 


8046. 


5 


2576. 


918.6 


5 


513.6 


5 


10240. 


4960. 


7 


1368. 


452.4 


7 


298.9 


7 


4119. 


1840. 


9 


838.8 


256.5 


9 


195.3 


9 


1875. 


820.2 


11 


545.2 


162.7 


11 


143.4 


11 


986.3 


413.9 


13 


371.4 


108.3 


13 


103.9 


13 


554.6 


232.4 


15 


273.5 


78.46 


15 


78.28 


15 


337.7 


137.8 


17 


206.6 


55.28 


17 


60.18 


17 


226.5 


88.37 


19 


162.1 


45.82 


19 


51.11 


19 


162. 


57.77 


21 


130.4 


33.90 


21 


40.63 


21 


107.4 


41.12 


23 


102.5 


26.90 


23 


34.76 


23 


71.90 


28.14 


25 


84.26 


22.64 


25 


28.13 


25 


59.46 


21.23 


27 


70.85 


18.27 


27 


24.60 


27 


38.62 


15.25 


29 


60.26 


15.58 


29 


21.12 


29 


30.19 


12.05 


31 


51.43 


13.08 


31 


17.05 


31 


25.45 


8.619 


33 


45.92 


11.02 


33 


17.66 


33 


22.84 


6.839 


35 


40.26 


9.718 


35 


16.21 


35 


15.55 


5.642 


37 


33.92 


7.860 


37 


12.86 


37 


13.24 


4.484 


39 


29.80 


7.281 


39 


10.61 


39 


11.64 


3.454 



Table 6: 
[The rapidity-integrated p\ distribution is given for single charm (NLO and Born) 
and the p\ distribution in the range \y\ < 1 is given for cc pair production (NLO 
only). The tabulated results have not been corrected for the rapidity bin width. 
The rapidity-integrated pair mass distribution is also given. All distributions are at 
sfs = 14 TeV and calculated with MRS D-' parton distributions. Note the 2 GeV bin 
width for the distributions.! 



bb Production y/s = 200 GeV 



dab/dp"^ 


(/ib/2 GeV 2 ) 


da bl /dp 2 T dy (/ib/2 GeV 2 ) 


da bl /dM (/ib/2 


GeV) 


V \ (GeV 2 ) 


NLO 


LO 


V \ (GeV 2 ) 


NLO 


M (GeV) 


NLO 


LO 


1 


0.2201 


0.1123 


1 


0.2073 








3 


0.1704 


0.0883 


3 


0.0524 








5 


0.1558 


0.0680 


5 


0.0263 








7 


0.1064 


0.0541 


7 


0.0170 








9 


0.1035 


0.0577 


9 


0.0118 


9 


0.0463 


0.0320 


11 


0.0863 


0.0406 


11 


0.00814 


11 


0.4363 


0.2100 


13 


0.0605 


0.0343 


13 


0.00660 


13 


0.3184 


0.1640 


15 


0.0478 


0.0255 


15 


0.00441 


15 


0.1987 


0.1050 


17 


0.0458 


0.0264 


17 


0.00341 


17 


0.1225 


0.0637 


19 


0.0351 


0.0190 


19 


0.00311 


19 


0.0753 


0.0400 


21 


0.0359 


0.0186 


21 


0.00274 


21 


0.0492 


0.0249 


23 


0.0300 


0.0139 


23 


0.00237 


23 


0.0318 


0.0160 


25 


0.0244 


0.0122 


25 


0.00201 


25 


0.0214 


0.0104 


27 


0.0216 


0.0116 


27 


0.00183 


27 


0.0145 


0.00688 


29 


0.0202 


0.0103 


29 


0.00156 


29 


0.0091 


0.00466 


31 


0.0171 


0.0080 


31 


0.00147 


31 


0.0069 


0.00321 


33 


0.0159 


0.0083 


33 


0.00121 


33 


0.0047 


0.00215 


35 


0.0125 


0.0054 


35 


0.00111 


35 


0.0032 


0.00154 


37 


0.0101 


0.0055 


37 


0.00111 


37 


0.0022 


0.00108 


39 


0.0097 


0.0049 


39 


0.00086 


39 


0.0016 


0.00075 



Table 7: 
[The rapidity-integrated p\ distribution is given for single b quarks(NLO and Born) 
and the °p\ distribution in the range \y\ < 0.35 is given for bb pair production (NLO 
only). The tabulated results have not been corrected for the rapidity bin width. 
The rapidity-integrated pair mass distribution is also given. All distributions are at 
y/s = 200 GeV and calculated with MRS D-' parton distributions. Note the 2 GeV 
bin width for the distributions.! 



bb Production y/s = 500 GeV 



dab/dp"^ 


(/ib/2 GeV 2 ) 


daJdp 2 T dy (/ib/2 GeV 2 ) 


daJdM (/ib/2 GeV) 


V \ (GeV 2 ) 


NLO 


LO 


V \ (GeV 2 ) 


NLO 


M (GeV) 


NLO 


LO 


1 


0.9809 


0.4798 


1 


0.3427 








3 


0.7911 


0.4024 


3 


0.2503 








5 


0.6490 


0.3362 


5 


0.1260 








7 


0.5492 


0.2801 


7 


0.0818 








9 


0.4528 


0.2358 


9 


0.0558 


9 


0.2652 


0.1199 


11 


0.3807 


0.1987 


11 


0.0426 


11 


1.737 


0.8711 


13 


0.3256 


0.1688 


13 


0.0341 


13 


1.436 


0.7552 


15 


0.2781 


0.1433 


15 


0.0285 


15 


0.9909 


0.5222 


17 


0.2428 


0.1248 


17 


0.0235 


17 


0.6646 


0.3503 


19 


0.2068 


0.1057 


19 


0.0197 


19 


0.4547 


0.2356 


21 


0.1824 


0.0932 


21 


0.0169 


21 


0.3132 


0.1612 


23 


0.1595 


0.0811 


23 


0.0147 


23 


0.2183 


0.1121 


25 


0.1429 


0.0719 


25 


0.0133 


25 


0.1566 


0.0797 


27 


0.1240 


0.0622 


27 


0.0122 


27 


0.1126 


0.0578 


29 


0.1108 


0.0557 


29 


0.0109 


29 


0.0850 


0.0419 


31 


0.0984 


0.0492 


31 


0.0098 


31 


0.0640 


0.0315 


33 


0.0898 


0.0435 


33 


0.0085 


33 


0.0469 


0.0236 


35 


0.0789 


0.0387 


35 


0.0076 


35 


0.0367 


0.0179 


37 


0.0716 


0.0350 


37 


0.0071 


37 


0.0291 


0.0138 


39 


0.0646 


0.0319 


39 


0.0074 


39 


0.0220 


0.0108 



Table 8: 
[The rapidity-integrated p\ distribution is given for single b quarks(NLO and Born) 
and the °p\ distribution in the range \y\ < 0.35 is given for bb pair production (NLO 
only). The tabulated results have not been corrected for the rapidity bin width. 
The rapidity-integrated pair mass distribution is also given. All distributions are at 
y/s = 500 GeV and calculated with MRS D-' parton distributions. Note the 2 GeV 
bin width for the distributions.! 



bb Production y/s = 5.5 TeV 



da b /dp 2 T (/ib/2 GeV 2 ) 


da bl /dp 2 T dy (/ib/2 GeV 2 ) 


da hl jdM (/ib/2 


GeV) 


V \ (GeV 2 ) 


NLO 


LO 


V \ (GeV 2 ) 


NLO 


M (GeV) 


NLO 


LO 


1 


23.59 


11.22 


1 


-2.366 








3 


19.38 


9.650 


3 


12.80 








5 


16.25 


8.253 


5 


6.634 








7 


13.84 


7.028 


7 


4.424 








9 


11.83 


6.065 


9 


3.303 


9 


6.102 


2.498 


11 


10.14 


5.148 


11 


2.496 


11 


42.57 


19.58 


13 


8.916 


4.469 


13 


1.946 


13 


37.41 


18.51 


15 


7.776 


3.890 


15 


1.726 


15 


27.66 


13.89 


17 


6.883 


3.424 


17 


1.439 


17 


20.00 


9.930 


19 


6.132 


3.004 


19 


1.199 


19 


14.41 


7.187 


21 


5.436 


2.650 


21 


1.073 


21 


10.53 


5.190 


23 


4.825 


2.296 


23 


0.9512 


23 


8.007 


3.863 


25 


4.357 


2.098 


25 


0.8151 


25 


6.028 


2.911 


27 


3.959 


1.875 


27 


0.7535 


27 


4.583 


2.202 


29 


3.545 


1.666 


29 


0.6718 


29 


3.577 


1.721 


31 


3.208 


1.526 


31 


0.5796 


31 


2.879 


1.342 


33 


2.950 


1.367 


33 


0.5276 


33 


2.248 


1.078 


35 


2.683 


1.207 


35 


0.5491 


35 


1.813 


0.8730 


37 


2.468 


1.131 


37 


0.4692 


37 


1.507 


0.7100 


39 


2.255 


1.034 


39 


0.4334 


39 


1.261 


0.5682 



Table 9: 
[The rapidity-integrated p\ distribution is given for single b quarks(NLO and Born) 
and the p\ distribution in the range \y\ < 1 is given for bb pair production (NLO 
only). The tabulated results have not been corrected for the rapidity bin width. 
The rapidity-integrated pair mass distribution is also given. All distributions are at 
y/s = 5.5 TeV and calculated with MRS D-' parton distributions. Note the 2 GeV 
bin width for the distributions.! 



bb Production y/s = 14 TeV 



da b /dp 2 T (/ib/2 GeV 2 ) 


da bl /dp 2 T dy (/ib/2 GeV 2 ) 


daJdM (/ib/2 GeV) 


V \ (GeV 2 ) 


NLO 


LO 


V \ (GeV 2 ) 


NLO 


M (GeV) 


NLO 


LO 


1 


68.43 


32.54 


1 


-13.36 








3 


56.73 


28.24 


3 


34.99 








5 


47.74 


24.25 


5 


17.94 








7 


41.32 


20.92 


7 


11.83 








9 


35.45 


18.10 


9 


8.519 


9 


17.57 


6.876 


11 


30.61 


15.55 


11 


6.833 


11 


124.0 


55.90 


13 


27.07 


13.60 


13 


5.537 


13 


112.4 


54.74 


15 


23.97 


11.93 


15 


4.665 


15 


85.11 


42.17 


17 


21.22 


10.41 


17 


3.813 


17 


62.92 


30.97 


19 


18.86 


9.192 


19 


3.392 


19 


46.41 


22.58 


21 


16.84 


8.225 


21 


3.125 


21 


34.27 


16.62 


23 


15.20 


7.227 


23 


2.618 


23 


26.12 


12.44 


25 


13.71 


6.477 


25 


2.328 


25 


19.89 


9.457 


27 


12.61 


5.878 


27 


2.112 


27 


15.51 


7.304 


29 


11.20 


5.215 


29 


1.772 


29 


11.93 


5.673 


31 


10.43 


4.710 


31 


1.811 


31 


9.610 


4.538 


33 


9.520 


4.368 


33 


1.588 


33 


7.908 


3.587 


35 


8.651 


3.962 


35 


1.409 


35 


6.267 


2.966 


37 


7.795 


3.492 


37 


1.349 


37 


5.132 


2.402 


39 


7.272 


3.245 


39 


1.279 


39 


4.323 


2.017 



Table 10: 
[The rapidity-integrated p\ distribution is given for single b quarks(NLO and Born) 
and the p\ distribution in the range \y\ < 1 is given for bb pair production (NLO 
only). The tabulated results have not been corrected for the rapidity bin width. 
The rapidity-integrated pair mass distribution is also given. All distributions are at 
yfs = 14 TeV and calculated with MRS D-' parton distributions. Note the 2 GeV bin 
width for the distributions.! 



Figure Captions 



1. (a) Gluon distributions from GRV HO (solid), MRS DO' (dashed), MRS D-' (dot- 
dashed) at Q = 2.4 GeV and GRV HO (dotted) at Q = 1.2 GeV. (b) The running of 
the coupling constant with scale. 

2. Investigation of uncertainties in the total cross section as a function of scale. 
Variation of the cc production cross sections with scale at (a) RHIC and (b) LHC. 
Variation of the bb production cross sections with scale at (c) RHIC and (d) LHC. 
Variation of the cc production cross sections at y/s at 200 GeV with jir at fixed [if (e) 
and with /xp at fixed lir (b). In each case, the circles represent the NLO calculation, 
the crosses, the Born calculation. 



3. Total charm production cross sections from pp and pA measurements [18], [19], |20 
pl| , |22| compared to calculations. The curves are: MRS D-' m c = 1.2 GeV, \x = 2m c 
(solid); MRS DO' m c = 1.2 GeV, \x = 2m c (dashed); GRV HO m c = 1.3 GeV, fi = m c 
(dot-dashed); GRV HO m c = 1.5 GeV, Li = m c (dotted). 



4. (a) Total charm production cross sections from n~p measurements jig , |25| , |2q , 
|27|, |8| compared to calculations. The curves are: GRV HO m c = 1.3 GeV, ii = m c 
on a nucleon (solid) and proton target (dashed); MRS D-' m c = 1.2 GeV, /x = 2m c 
with SMRS P2 (dot-dashed) and SMRS PI (dotted) on a proton target, (b) The bb 
production cross section from n~p interactions pj| ^0], ^T| ^2|. The calculations use 
m b = 4.75 GeV and \i = m b . The curves use GRV HO (solid) and MRS D-' with 
SMRS PI (dashed). 



5. Comparison with D meson (a) p\ and (b) xf distributions at 800 GeV |T!|. The 



NLO calculations are with MRS D-' (solid) and GRV HO (dashed) parton distribu- 
tions. 

6. Comparison with DD production for (a) p\ and (b) M and (c) <\> at 800 GeV 
[ J36fl . The NLO calculations are with MRS D-' (solid) and GRV HO (dashed) parton 
distributions. 



7. Comparison with b quark production cross sections at (a) UA1 |37| and (b) CDF 
[ J3S|| . The NLO calculations are with MRS D-' (solid) and GRV HO (dashed) parton 
distributions. 

8. Predictions for c and cc production at y/s = 200 GeV with MRS D-' distributions. 
The c quark px distributions at NLO (solid) are shown in (a) and the rapidity distri- 
butions at LO (dashed) and NLO (solid) are shown in (b). The cc pair distributions 
are shown in (c)-(f). The LO (dashed) distributions are shown only for mass and 
rapidity. Additionally, the px and px distributions are shown with a central cut in 
rapidity. (The rapidity bin widths are removed.) The corresponding distributions 



from HIJING are also shown, again with the rapidity bin width divided out. 

9. Predictions for c and cc production at y/s = 200 GeV with GRV HO distributions. 
The corresponding distributions from HIJING are also shown. 

10. Predictions for c and cc production at y/s = 500 GeV with MRS D-' distributions. 

11. Predictions for c and cc production at y/s = 500 GeV with GRV HO distributions. 

12. Predictions for c and cc production at y/s = 5.5 TeV with MRS D-' distributions. 
The corresponding distributions from HIJING are also shown. 

13. Predictions for c and cc production at y/s = 5.5 TeV with GRV HO distributions. 
The corresponding distributions from HIJING are also shown. 

14. Predictions for c and cc production at y/s = 14 TeV with MRS D-' distributions. 

15. Predictions for c and cc production at y/s = 14 TeV with GRV HO distributions. 

16. Predictions for b and bb production at y/s = 200 GeV with MRS D-' distributions. 

17. Predictions for b and bb production at y/s = 200 GeV with GRV HO distributions. 

18. Predictions for b and bb production at y/s = 500 GeV with MRS D-' distributions. 

19. Predictions for b and bb production at y/s = 500 GeV with GRV HO distributions. 

20. Predictions for b and bb production at y/s = 5.5 TeV with MRS D-' distributions. 

21. Predictions for b and bb production at y/s = 5.5 TeV with GRV HO distributions. 

22. Predictions for b and bb production at y/s = 14 TeV with MRS D-' distributions. 

23. Predictions for b and bb production at y/s = 14 TeV with GRV HO distributions. 



24. Dilepton (a) mass, (b) rapidity, and (c) p? distributions at ^/s = 200 GeV from 
cc decays calculated using MRS D-' distributions are shown. 

25. Dilepton distributions at ^/s = 200 GeV from bb decays. 

26. Dilepton distributions at y/s = 5.5 TeV from correlated cc decays. 

27. Dilepton distributions at \/s = 5.5 TeV from uncorrelated cc decays. 

28. Dilepton distributions at y/s = 5.5 TeV from bb decays. 
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